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ABSTRACT 

Suppose an astronomer is equipped with a device capable of detecting emissions — 
whether they be electromagnetic, gravitational, or neutrino — from transient sources 
distributed throughout the cosmos. Because of source rate density evolution and vari- 
ation of cosmological volume elements, the sources first detected when the machine is 
switched on are likely to be ones in the high-redshift universe; as observation time in- 
creases, rarer, more local, events will be found. We characterize the observer's evolving 
record of events in terms of a 'probability event horizon', converging on the observer 
from great distances at enormous speed, and illustrate it by simulating neutron star 
birth events distributed throughout the cosmos. As an initial application of the con- 
cept, we determine the approach of this horizon for gamma-ray bursts (GRBs) by 
fitting to redshift data. The event rates required to fit the model are consistent with 
the proposed link between core-collapse supernovae and a largely undetected popula- 
tion of faint GRBs. 
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1 INTRODUCTION 

During the last decade, space-based and terrestrial astron- 
omy have revolutionized our understanding of the high-z 
Universe (redshift z > 1). Astronomers can now observe 
the high-z Universe in many bands of the electromagnetic 
spectrum and have begun putting upper limits on the gravi- 
tational wave spectrum. Their telescopes and detectors have 
provided a direct probe of cataclysmic transient astrophysi- 
cal phenomena that are rare in our local Universe but occur 
at high rates on a cosmological scale. For example, optical 
observations of Type la supernovae (SNe), are used as tools 
to constrain the cosmological model that best describes our 
Universe. Their uniformity in luminosity means they can 
be used as standard candles for determining cosmological 
distances. Recently 23 high-z Type la SNe, including 15 at 
> 0.7 were used as a cosmic ruler to show that the Universe 
is accelerating and h as cosmological pa rameters consistent 
with spatial flatness (Barris, 

ieTaLllio^)- 

The high-energy signature of one class of event occurs in 
the gamma-ray region of the electromagnetic spectrum: Cos- 
mological gamma-ray bursts (GRBs) are extremely energetic 
transient events and their fluence (time-integrated flux) and 
time variability show that their power sources must be com- 
pact. At least a fraction of GRBs are associated with massive 
and short-lived progenitor stars, implying that either a neu- 
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tron star or black hole must be driving the emissions. Recent 
estimates suggest that the universal cosmological GRB rate 
could be as high as one per minute ( van Putton & Rcgimb^ 
l2004h. 



We adopt these deflnitions: an event is a cataclysmic 
astrophysical transient occurrence, such as a core-collapse 
SN or a GRB, with a duration much less than the observa- 
tion period; an observer records event times and redshifts, 
obtained as data from a flux-limited detector sensitive to 
high-z events. 

With the launch of the Swift satellite (on 2004 
November 20), a multi- wavelength GRB observatory 
(http:/ /swift. gsfc.nasa.gov/ ), up to 1000 GRBs are expected 
to be detected in 3 years as a result of a sensitivity limit 5 
times fainter than that of BATSE (the Burst And Transient 
Source Experiment that was aboard the Compton Gamma 
Ray Observatory). The Swift mission has the potential to 
localize about one GRB per day, providing the possibility of 
cataloging hundreds of GRB redshifts. 

Predicting the rates of transient astrophysical events, 
such as GRBs, is fundamental to detector design and plan- 
ning realistic science goals. An understanding of the proba- 
bility of achieving those goals is correlated with the proba- 
bility of acquiring a useful sample of events. The 'ultimate 
detector' would be sensitive enough to detect all events out 
to the redshift where the events of the relevent type flrst 
started. For example, determining the event rate for GRBs 
as a function of z would ideally require detectors with all-sky 
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coverage and extraordinary sensitivity. In addition, selection 
effects plague many high-z observations. 

With continuing technical improvements to astronomi- 
cal instrumentation, it is conceivable that the observed rates 
of GRBs will approach the theoretical limit in the future. 
But for the present, all astronomical detectors probe only a 
fraction of the total volume encompassing all possible events. 
For any astronomical detector and source type, one can de- 
fine a 'detectability horizon' centred on the detector and 
encompassing the volume in which such events are poten- 
tially detectable; the horizon distance is determined by the 
flux limit of the detector and the source flux. A second hori- 
zon, defined by the minimum distance for at least one event 
to occur over some observation time, with probability above 
some selected threshold, can also be defined. We call this 
the 'probability event horizon' (PEH). It describes how an 
observer and all potentially detectable cosmological events 
of a particular type are related via a probability event dis- 
tribution encompassing all such events. 



2 PROBABILITY EVENT HORIZON 

2.1 PEH in a Euclidean universe 

To illustrate the idea, we first show how a PEH can be de- 
fined in the simple case of a static Euclidean universe, as- 
suming an isotropic and homogenous distribution of events. 
For a constant event rate ro per unit volume, the mean 
cumulative event rate in a volume of radius r is R{r) = 
(4/3)7rr^ro. The events are independent of each other, so 
their distribution is a Poisson process in time: the proba- 
bility for at least one event to occur in this volume during 
observation time T at a mean rate R{r) at constant proba- 
bility e is given by an exponential distribution: 

p(n>l;J?(r),r) = l-e'^M^ = e, (1) 

e~^"^ being the probability of zero events occurring. For Eq. 

to remain satisfied as T increases, the mean number of 
events in this volume, A^e ~ R{r)T = |ln(l — e)|, must remain 
constant. The corresponding radial distance, a decreasing 
function of observation time, defines the PEH: 



(2) 



The speed at which the horizon approaches the observer — 
the PEH velocity — is obtained by differentiating (T) with 
respect to T: 



(3) 



In practice, we'll take e = 0.95, corresponding to a 95% 
probability of observing at least one event within distance r, 
and so A^e = 3 is the mean number of events. Note that rf 
and wf^^ depend on the choice of e through N^''^: if e = 0.1 
is chosen, =0.1 and rf^^ and vf^^ will be bigger by a 
factor of 3 than for e = 0.95. 

If we consider event rates in volumes of cosmological 
scale, vf^^{T) can be very much greater than the speed 
of light. For neutron star formation, ro « 5 x lO"^'^ s~^ 
Mpc~^ in the nearby Universe (Coward, Burman & Blair 
2001, section 2.1). Taking this value and solving Eq. for 
T with vf^^{T) =c = 9.69 x lO"^'^ Mpc s"^ and e = 0.95, 
we find that about 2 x 10^ years of observation time elapse 



before vf^^{T) drops below c. The horizon would then be 
about 0.2 Mpc from the observer, but the assumption of 
homogeneity breaks down for distances less than a few tens 
of Mpc. After one year of observation, the PEH would be 30 
Mpc from the observer, which is about the distance within 
which we would expect fluctuations in R{r) resulting from 
the local distribution of galaxies to dominate; wf^^(r) at 
this distance is about 2 pc s~^. 



2.2 Cosmological factors 

To extend the PEH model to a Friedmann cosmology re- 
quires modifying the simple Euclidean model to incorporate 
cosmic evolution. The expansion of the Universe is described 
by the evolving Hubble parameter H{z), which can be ex- 
pressed in terms of th e contributions of matter and vacuum 
energy (|Peeblejll99l : 



h{z) = H{z)/Ho = [nn.{l + zf + ^a] 



1/2 



(4) 



for a 'flat-A' cosmology (a spatially flat cosmology with cos- 
mological constant). 

The comoving radial (line-of-sight) distance is obtained 
from h{z) by 



(5) 



this corresponds to integrating elements of comoving 
distance — ones that remain constant with epoch for neigh- 
bouring objects moving with the Hubble flow — and is the 
appropriate distance me asure for structures locked to the 
Hubble flow l! Hogdl2000l) . 

The luminosity distance, Dl: deflned so that the 
inverse- square law of intensity applies, it can be measured 
by means of 'standard candles'. These two distance measures 
are related through redshift by Dl = (1 + z)Dc. 

Volume elements in a Friedmann cosmology can be ex- 
pressed in terms of Dc{z) and H{z) (e.g. Porciani & Madau 
2001): 



dz Ho h(z) 



(6) 



In a flat-A cosmology, the density parameters of mat- 
ter and vacuum energy sum to unity; we use fim = 0.3 
and r^A = 0.7 for their present-epoch values and take 
Ho = 70 km s~^ Mpc~^ for the Hubble parameter a,t z = 
0. For comparison, we also consider the Einstein-de Sit- 
ter (EdS) cosmology — a flat Universe with zero vacuum 
energy — in which f2m = 1 and Qa ~ 0. Figure^plots dV/dz 
as a function of comoving radial distance. It shows that vol- 
ume elements (of fixed Az) for the EdS and flat-A cosmolo- 
gies reach a maximum and then decline, a result of cosmic 
expansion. 



2.3 Source rate density evolution 

Neutron star (NS) births result from massive short-lived pro- 
genitor stars and there is mounting evidence that GRBs do 
so too, so both types of events should closely track the evolv- 
ing star formation rate (SFR). 

We use a dimensionless SFR density evolution factor 
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Figure 1. The dimensionless volume element per unit redsiiift, 
dV/dz (cf Eq. (6) with c = 1 = Hq) as a function of comoving 
radial distance for the Euclidean, flat-A (0.3, 0.7), and EdS cos- 
mologies; h{z) = 1 for the Euclidean case. Here we use shells of 
distance rather than shells of z because redshift is inapplicable 
for the static Euclidean cosmology. The curves show that event 
counts in the flat-A and EdS cosmologies reach a maximum and 
then decline. This is modified by source rate density evolution 
and, for event rates, by time dilation from cosmic expansion. For 
distances less than about 1 Gpc the curves converge to that of 
the Euclidean cosmology. 



e{z), normalized to unity in our local intergalactic neigh- 
borhood, to account for source rate evolution: 

(1 + M^)e'3" 



+ W 



(7) 



The two SFR models used in this paper have the parameter 
values 

(i) SFl: {Q,R,W) = (3.4,3.8,45) 

(ii) SF2: {Q,R,W) = (3.4,3.4,22) 

as fitted u sing the EdS cosmology vyith Ho = 65 km s^^ 
Mpc"^ bv lPorciani fc Madaul 1I2OOIII . The first star forma- 
tion history matches most Ha luminosity densities measured 
in the UV, and includes an upward correction for dust red- 
dening. The second model includes more substantial dust 
extinction at high 2. 

An observation-based SFR density is cosmology- 
dependent, because it is deduced from the observed spec- 
tral luminosity in set volumes. The above SFR models can 
be re-scaled to other cosmolo gies, as shown in the appendix 
of IPorciani fc Madaul i200ll) . The luminosity density, and 
hence SFR density, are proportional to Di^'^ / (dV / dz) , which 
is proportional to (1 -I- z)^H{z), and so both scale as H{z) 
among different cosmologies. But e{z), which is normalized 
to unity ai z = 0, just follows the shape of the SFR density 
curve and so scales as h{z); that is e{z)/h{z) is invariant un- 
der change of cosmology. For converting the above evolution 
factors to a flat-A cosmology with Qm ~ 0.3 and CIa = 0.7, 
the scale factor is 



[0.3{1 + z f + 0.7]^^^ {1 + z) 



-3/2 



(8) 



as follows from Eq. Q for h{z). The resulting e{z) functions 
are shown in Fig.|5| 
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Figure 2. Dimensionless evolution factors, e(z), normalized to 
unity at 2 = 0, based on the parame trized star formation rate 
scenarios labelled as SFl and SF2 in iPorciani fc Madaul J200ll) 
but scaled to the fiat-A (0.3, 0.7) cosmology; model SF2 allows 
for high dust extinction at high z. 
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Figure 3. The differential neutron star birth rate, d-R/dz, as 
a function of redshift, for the flat-A (0.3, 0.7) and EdS cos- 
mologies and the star formation rate model labelled SFl in 
[Porciani fc Madaul ^2001^ but scaled to the flat-A (0.3, 0.7) cos- 
mology. The present-epoch rate ro is taken to be 5 X 10~^^ s~^ 
Mpc-3. 



2.4 PEH in flat-A and EdS cosmologies 

In the Euclidean model, the event rate in the shell r to 
r dr centred on the observer is dR = 47rr^rodr. In the 
standard Friedmann cosmologies, one can express the differ- 
ential event rate as the event rate in the redshift shell z to 
z -I- dz: 



diJ^^^^d., 

dz \ + z 



(9) 



where dV is the cosmology-dependent co-moving volume el- 
ement and R(z) is the all-sky (47r solid angle) event rate, 
as observed in our local frame, for sources out to redshift z. 
Source rate density evolution is accounted for by the dimen- 
sionless evolution factor e(z)\ as this is normalized to unity 
in the present-epoch universe (2; = 0), ro is now the 2 = 
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Figure 4. The cumulative neutron star birth rate, R{z), obtained 
by integrating dR/dz from 2 = 0. Universal event rates (inte- 
grated throughout the cosmos) are 14 and 27 in the EdS and 
flat-A cosmologies respectively. 

rate density. The (1 + z) factor accounts for the time dila- 
tion of the observed rate by cosmic expansion, converting a 
source-count equation to an event rate equation. 

The cumulative NS birth or GRB event rate R[z) is cal- 
culated by integrating the differential rate from the present 
epoch to redshift z. It depends on cosmology through the 
factors AV/Az and e{z) in Eq. for d_R. The differential 
and cumulative rates of NS formation as functions of z are 
plotted in Figs.l^land^] for model SFl. Figure^Jshows that 
R{z) is higher in the flat-A (0.3, 0.7) cosmology than in the 
EdS one. Because the curves scale with ro, they can also be 
applied to GRBs, but with the rates a factor of about 10^ 
lower for the classical GRBs. 

For event types locked to the star formation rate, e{z) 
is obtained by normalizing a SFR model to its z = value. 
Presently there is no consensus among astronomers on a 
model of star formation history at high z, so we shall use 
either a c onstant star formation hist ory or the parametrized 
models of lPorciani fc Madaul (1200 ll) . 

Integrating the differential rate yields the cosmological 
analogue of the Euclidean mean cumulative event rate: 

R{z) = [ {AR/Az) Az. (10) 
Jo 

Because of variation of cosmological volume elements 
(Fig.l), the rate derivative AR/Az reaches a maximum at 
high z in both the flat-A and EdS cosmologies. Consequently 
there is a slower growth in the cumulative rate as compared 
to that in the simple Euclidean model in which the volume 
elements grow as . 

The events follow the same probability distribution as 
in the Euclidean case (Eq. but with R now expressed as 
a function of z: 

p(n> l;i?(z),r) = l-e-"(")^ = e, (11) 

with Nt = R{z)T — jln(l — e)j constant. Instead of inverting 
to find rf^"(r), we flnd zf^^{T) by solving this condition 
numerically, thus defining the PEH. Converting zf^^{T) to 
luminosity distance and difi^erentiating with respect to T 
then yields ^^'^"(r). 




Observer time (s) 



Figure 5. The probability event horizon (PEH) velocity for Eu- 
clidean, flat-A (0.3, 0.7) and EdS cosmologies assuming a constant 
event rate density of 5 X lO"^'^ s~^ Mpc"''. The starting points at 
the upper left indicate that the PEH is 'stalled' until the number 
of events exceeds some threshold, which is cosmology dependent 
as a result of different cumulative event rate functions. The flat-A 
and EdS curves depart markedly from the Euclidean curve dur- 
ing the first seconds of observation time. The change in gradient 
apparent in the flat-A and EdS curves occurs because of a slower 
increase in volume at high z. 




Redshift, z 

Figure 6. Contours, based on the same constant event rate den- 
sity and flat-A cosmology as used in Fig. |5] show the probability 
(labelling the curves) for at least one event to occur in the redshift 
shell z to z + dz a.s a, function of z and observation time. The min- 
imum of each of the contours corresponds to the most probable 
redshift at which an event will occur. As a function of obser- 
vation time, each probability contour follows two trajectories — 
one approaching the observer and the other receding. The reced- 
ing contours reach the redshift where events first started within 
a fraction of a second, as opposed to years in observer time for 
an approaching probability contour to reach the local Universe. 

3 NEUTRON STAR BIRTH EVENTS 

Figure shows the PEH velocity for NS birth events in the 
three cosmologies — Euclidean, fiat-A and EdS — with no 
source evolution so e{z) — 1. Because the local Universe is 
approximately Euclidean, the PEH velocity converges to the 
Euclidean result when the observation time becomes much 
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Figure 7. Simulation illustrating the rapid motion of the PEH 
for the first few tens of seconds of observation time, and a model 
curve with a probability threshold e = 0.95 for at least one event 
to occur at a distance less than that of the PEH. We assume 
a universal cumulative event rate of about 25 s"'^ as seen in our 
frame, comparable to the NS birth rate integrated throughout the 
Universe, a fiat-A (0.3, 0.7) cosmology and a source rate evolution 
locked to the SFR model labelled SF2 in Porciani & Madau (2001) 



longer than the mean time interval between events. The cu- 
mulative rate, R{z), is cosmology dependent, because the 
differential rate is determined by the evolution of volume 
elements as a function of z (Eq. |SJ . Variations in event rates 
result in differing PEH velocities; for example, the relatively 
smaller event rate obtained using the EdS model (see Fig. 
|1J manifests as a longer time delay before the initiation of 
PEH motion because p{n ^ 1; R{z),T) remains below 0.95 
for a longer observation time. 

Figure |S| shows probability contours, in the flat- A cos- 
mology, for at least one event to occur in the redshift shell 
z to z + dz as a. function of T and z. The curves reveal two 
horizons, one approaching and one receding as T increases; 
the receding horizon, corresponding to the observation of in- 
creasingly distant events beyond the peak of d_R/dz, is very 
short lived. 

To demonstrate the motion of a PEH we have simu- 
lated the occurrence, as a function of observation time, of 
NS birth events distributed throughout the Universe. The 
redshifts are treated as a random variable following a prob- 
ability distribution function obtained by normalizing dR/dz 
(Coward, Burman & Blair 2002a,b): 



p{z)Az = AR/R{oo) 



(12) 



R{oo) being the Universal rate, integrated throughout the 
cosmos, as seen in our frame. The corresponding cumulative 
distribution function C{z), giving the probability of an event 
occurring in the redshift range to z, is 



C{z) = R{z)/R{oo) 



(13) 



We use the inverse function of C{z) to produce values of z 
by employing a random number generator to select values 
from C{z). 

Our simulation uses ro ~ 5 X 10"^^ s"^ Mpc"^ for the 
z — rate density of these events and an evolution locked 
to the SFR model labelled SF2 in Porciani & Madau (2001) 
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Figure 8. Filled diamonds show the closest events as a function 
of observation time for the 34 GRBs for which the redshift is 
known with reasonable confidence; the solid square is the highly 
under-luminous GRB 980425. The observations span a 7-yr pe- 
riod, starting from the first redshift measurement in 1997 Febru- 
ary and continuing to 2004 July. The solid curve shows the PEH, 
excluding the outlier GRB980425, assuming a flat-A (0.3, 0.7) cos- 
mology and a source rate evolution based on the star formation 
rate model labelled SFl in Porciani &; Madau (2001). A local rate 
density ro ~ 0.8 yr~^ Gpc~^ is used to fit to the redshift data 
for the classical GRBs for the first 4-5 yr of observation. (The 
probability of GRB 980425 occurring inside a volume bounded 
by its redshift z = 0.0085 in this time is 5 X 10~^.) After this 
period, a more rapidly decreasing trend in distance is observed. 
The dashed curve shows the PEH model fitted to the data from 
2001 September to 2004 July with ro rti 27 yr"! Gpc^^*. GRBs 
associated with core-collapse SNe are labelled. 



scaled to the flat-A (0.3, 0.7) cosmology. The running mini- 
mum redshifts of the simulated events as a function of ob- 
servation time define the PEH; these have been converted 
to luminosity distances and plotted in Fig. 7 along with a 
PEH curve for e = 0.95. Of these 103 events only one is at 
a distance greater than the PEH. The binomial 95% confi- 
dence interval, give 0.9471 < P < 0.9998 for one failure in 
103 trials, which is in good agreement with e = 0.95 for at 
least one event to occur within the PEH. 

In this calculation of the approach of the PEH for sim- 
ulated NS birth events throughout the Universe, we have 
assumed that the observation time for each event is small 
compared with the typical observed temporal separation of 
events, which averages about 40 ms. The calculation is di- 
rected toward future 'Advanced LIGO'-type gravitational 
wave detectors, which may be able to detect a stochastic 
background of gravitational radiation generated by these 
events. Recent calculations, based on the computed wave 



functions of iDimmelmeier. Font fc MiilleJ (120021) . suggest 
that the duty cycle of the signal will be well below one (How- 
ell et al., 2004, figure 16), implying a signal consisting largely 
of non-overlapping wave-trains, described as 'popcorn noise'. 

For optical surveys of very high-z SNe, selection effects 
and the sensitivity limits of the detectors reduce the ob- 
served rate substantially. If these effects were well enough 
understood, then the rate density of events could still be in- 
ferred from the evolving record of observations. These ideas 
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Table 1. GRB rcdshifts adapted from the table of localized GRBs at 
liittp: //www. aip.de/~jcg/grbgen. html The instruments used to localize the GRBs are listed along 
with confirmation if the Interplanetary Network (IPN) was used. Other columns show: XA = X-ray 
afterglow; OT = optical (or IR) transient; RA = radio afterglow; lAUC (International Astronomical Union 
Circular) number; z = redshift. An 'X' after the date marks events with no gamma-ray detection, but with 
X-ray properties similar to GRBs — so-called X-ray flashes. The 'C' after one of the dates identifies it as 
the 3rd GRB on that day. The running minimum observed GRB redshifts from 1997 February to 2004 July 
(entries highlighted in bold) are used in Fig. El with GRB 980425 regarded as an outlier. 



GRB Instrument IPN XA OT RA lAUC 



970228 


SAX/WFC 


y 


y 


y 


n 


6572 


0.695 


970508 


bAX/WlC 




y 


y 


y 


6649, 6654 


0.835 


970828 


SAX/WFC 


y 


y 


n 


y 


6726, 6728 


0.9578 


971214 


SAX/WFC 


y 


y 


y 


n 


6787, 6789 


3.42 


980425 


SAX/WFC 




y 


SN 


y 


6884 


0.0085 


980613 


SAX/WFC 




y 


y 


n 


6938 


1.096 


980703 


RXTE/ASM 




y 


y 


y 


6966 


0.966 


990123 


SAX/WFC 


y 


y 


y 


y 


7095 


1.6 


990506 


BAT/PGA 


y 


y 




y 




1.3 


990510 


SAX/WFC 


y 


y 


y 


y 




1.619 


990705 


SAX/WFC 


y 


y 


y 


n 


7218 


0.86 


990712 


SAX/WFC 






y 


n 


7221 


0.434 


991208 


Uly/KO/NE 


y 




y 


y 




0.706 


000131 


Uly/KO/NE 


y 




y 






4.5 


000210 


SAX/WFC 


y 


y 




y 




0.846 


000301C 


ASM/Uly 


y 




y 


y 




2.03 


000911 


Uly/KO/NE 


y 




y 


y 




1.058 


000926 


Uly/KO/NE 


y 


y 


y 


y 




2.066 


010921 


HE/Uly/SAX 


y 




y 






0.45 


011121 


SAX/WFC 


y 


y 


y 


y 




0.36 


011211 


SAX/WFC 




y 


y 






2.14 


020124 


HETE 


y 




y 






3.198 


020405 


Uly/MO/SAX 


y 




y 


y 




0.69 


020813 


HETE 


y 


y 


y 


y 




1.25 


020903X 


HETE 






y 


y 




0.25 


021004 


HETE 




y 


y 


y 




2.3 


021211 


HETE 






y 






1.01 


030226 


HETE 




y 


y 






1.98 


030323 


HETE 








y 




3.372 


030328 


HETE 




y 


y 






1.52 


030329 


HETE 




y 


y 


y 


8101 


0.168 


030429 


HETE 






y 






2.65 


031203 


INTEGRAL 




y 


y 


y 


8250, 8308 


0.105 


040701X 


HETE 




y 








0.2146 



tion effects and the local rate density; these parameters are 
all uncertain. 

We use GRB redshift data from 1997 Febuary (when the 
first GRB redshift was obtained) to 2004 July, which consist 
of 235 GRBs localized to less than 1 degree within a few 
hours to days. About 40 GRB redshifts have been measured 
but only 34 are known with certainty (mostly identified from 
their host galaxies); these events, along with details of each 
observation, are shown in Table 1. 

Figure |5] plots the temporal evolution of the minimum 
observed luminosity distance, using these measured GRB 
redshifts, as a function of observation time from 1997 Febru- 
ary to 2004 July; we do not include GRB redshifts that are 
highly uncertain. We note that the very under-luminous and 
nearby GRB 980425 (z = 0.0085 corresponding to 40 Mpc), 
and GRBs 030329 (2 = 0.168) and 031203 {z = 0.105), are 
associated with cc SNe, probably Type Ib/c. 

Our data points are obtained by selecting the sources 



will be demonstrated in the next section for GRB afterglow 
observations. 



4 GRB EVENTS 

The PEH velocity is based on the act of observation, where 
the observer is at the centre of a Universe defined in terms 
of an event probability distribution. The temporal evolution 
of the PEH for a particular class of cosmological transient 
event is defined by the history of event rates. Different types 
of events, for example NS births and GRBs, will have dif- 
ferent PEH velocities because their mean event rate densi- 
ties differ, even if both are locked to the SFR. The PEH 
concept can be applied to the observed GRB redshift distri- 
bution assuming some knowledge of source-rate evolution, 
luminosity distribution of the sources, observational selec- 
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with minimum redshift as a function of observation time 
over the GRB redshift data set. A PEH model curve is fit- 
ted to these data, assuming the flat-A (0.3, 0.7) cosmology 
and a s ource rate evolution locked to the SFR model labelled 
SFl in IPorciani fc Madaul (|20o3). The differential rate has 
one free parameter, the present-epoch GRB rate density, ro, 
which we vary in the PEH model to fit the data. We scale the 
differential event rate by a GRB luminosity function follow- 
ing a log-normal distribution (jBromm fc Loeb 2002), with a 
mean GRB luminosity of 2 x 10^® photons s"'^, and assume 
a photon flux threshold of 0.2 s"^ cm~^ for BATSE (see Ap- 
pendix A). We note that the form of luminosity distribution 
that would fit both the 'classical' GRBs and GRBs 980425 
and 031203 is not apparent. 

Excluding the outher GRB 980425, we find that the 
data for the first 4-5 yr of observation time (classical GRBs) 
are best fitted using a PEH curve with ro ~ 6 x 10^^ 
yr~^ Gpc""^. Scaling this rate by accounting for the solid 
angle monitored by t he BeppoSAX camera or HETE tele- 
scope, about 0.9 str jGuettal liooil . yields ro ^ 0.8 yr"^ 
Gpc~^, similar to the rate density of 0.5 yr~^ Gpc~^ de- 
duced from observations of both optica l and optically dark 
GRBs llSchmiddl200ll iFrafl et alJl200ll) . The probability of 
at least one GRB occurring inside a volume bounded by 
z = 0.17 (in 4.5 yrs observation time) has the quite signifi- 
cant value of 0.37, but the probability of one occurring inside 
a volume bounded hy z — 0.0085 is only about 5 x 10~^: it 
is difficult to reconcile the very nearby G RB 980425 with 
the re dshift distribution of classical GRBs. ISoderberg et alJ 
Jio^ conclude that the two under- luminous GRBs 980425 
and 031203 were intrinsically sub-energetic and un-jetted 
events. 

Data from 2001-2004, which include GRBs 030329 and 
031203, show a more sharply decreasing trend in Fig. |S1 
implying a relatively higher PEH velocity and so a higher 
R{z) and hence a higher ro; fitting a PEH curve to this 
subset of data requires ro ~ 27 yr~^ Gpc~^. With this ro, 
the probability for an event to occur in the volume bounded 
by GRB 980425 {z = 0.0085) in 4.5 yr is 0.0016. Thus, 
it seems likely that GRB 980425 is a member of a more 
numerous population of low-luminosity GRBs. 

Because GRBs 980425 and 031203 were sub-energetic 
and showed no evidence for jetted emissions, ro for such a 
low-luminosity GRB population could be orders of magni- 
tude higher than that of the classical GRB popula tion, but 
still o nly a small fraction of the Type Ib/c SN rate iCowardI 
l2005l) . 



5 CONCLUDING REMARKS 

The probability event horizon concept provides a means of 
analyzing an observer's evolving record of events in terms 
of the distribution of the sources throughout the cosmos. It 
is particularly sensitive to the shape of the low probability 
'tail' of the source distribution. The dependence of the PEH 
velocity on the local source rate density means that by ana- 
lyzing the redshift distribution of the sources as a function of 
observation time one can infer the local source rate density. 

As an initial application of this, we have shown how 
fitting a PEH model to GRB redshift data can yield infor- 
mation on the local GRB rate density. Such an analysis can 



be used to probe the progenitor populations that comprise 
the distributions. This is demonstrated by Figure |H| where 
GRB 980425 is either a statistical outlier of the classical 
population or a member of a different distribution of GRBs 
that could be intrinsically faint; alternatively, it could be a 
classical GRB observed away from its jet axis. 

The PEH concept has applications to gravitational wave 
detection. Understanding the detectability of astrophysical 
sources as a function observation time is crucial to these 
searches. Current detectors are only sensitive to local events, 
such as binary coalescences out to several tens of Mpc; the 
PEH for these events corresponds to long observation times 
of the order of a year or so. Because of the rarity of such 
events, we are modelling the detectability of more frequent 
but fainter gravitational wave sources using the PEH con- 
cept. 

We are currently studying the application of the PEH 
in simulations of the evolving spectra of contributions to the 
stochastic gravitational wave background generated by NS 
births and other cataclysmic events distributed throughout 
the Universe. 
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APPENDIX A: GRB LUMINOSITY FUNCTION 

The GRB luminosity function (LF), together with the flux 
sensitivity threshold of the instrument, determines the frac- 
tion of a ll GRBs potentially detectable with that instrument 
jNorrij 12002 1. Following Bromm & Loob (2002), we scale 
AR/Az to account for the distribution of GRB luminosities: 



'(/'GRB (2)= / P{L)AL, (Al) 

where 'i/'GRB(z) is the GRB rate scaling function and p{Ij) 
is the GRB LF with L the intrinsic luminosity in units of 
photons s~^. With /um denoting the flux sensitivity thresh- 
old, in photons s~^ m~^, the minimum detectable luminos- 
ity can be expressed as a function of redshift by Lun^{z) = 
47r_DL^(2)/iim, with Di,{z) the luminosity distance. 

In modelling the GRB redshift dist r ibutio n we use the 
log- normal LF cited bv lBromm fc Loebl (|22q2)- 

where o and Lq are the width and average luminosity, re- 
spectively. We assume a = 2 and Lq — 2 x 10^^ s~^, with 
fiim = 0.2 photons s"^ cm"^ for BATSE jBromm fc Loebl 
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